Hydrodynamical instabilities in an expanding quark gluon plasma 
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We study the mechanism responsible for the onset of instabilities in a chiral phase transition at 
nonzero temperature and baryon chemical potential. As a low-energy effective model, we consider an 
expanding relativistic plasma of quarks coupled to a chiral field, and obtain a phenomenological chiral 
hydrodynamics from a variational principle. Studying the dispersion relation for small fluctuations 
around equilibrium, we identify the role played by chiral waves and pressure waves in the generation 
of instabilities. We show that pressure modes become unstable earlier than chiral modes. 
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I. INTRODUCTION 

It is commonly believed that QCD at very high tem- 
peratures and/or densities allows for a new phase of 
strongly interacting matter, the quark-gluon plasma 
(QGP). Compelling lattice QCD results point in that 
direction [lj, and experiments in ultra-relativistic heavy- 
ion collisions [2| under way at BNL-RHIC 3] and planned 
for CERN-LHC attempt to glimpse at this elusive state 
of matter that was presumably present in the early uni- 
verse 0. 

Depending on the nature of the QCD phase transition, 
the hadronization process of the expanding QGP gener- 
ated in a high-energy heavy-ion collision may proceed in 
a number of different ways |j|. In the analysis of phase 
conversion of strongly interacting matter from a chirally 
symmetric to a chirally broken state in the phase dia- 
gram of QCD, a mechanism of supercooling is usually 
implied [jj E|- The kinetics of domain formation and 
growth can then be described by the the mechanisms of 
nucleation or spinodal decomposition, depending on the 
degree of supercooling |{|. An expanding plasma will 
probe, in principle, different temperatures and experi- 
ment both mechanisms. In the case of the early universe, 
it is well known that the time scales for the expansion 
during the QCD transition are so slow that phase conver- 
sion is driven by the nucleation of bubbles of true vacuum 
inside the metastable phase, as the universe reaches tem- 
peratures below a critical value On the other hand, 
the QGP presumably created in experiments in ultra- 
relativistic heavy-ion collisions expands at a pace several 
orders of magnitude faster than the primordial universe 
H El , and might simply bypass nucleation, entering the 
domain of spinodal decomposition H H [H El El • 
Some results from CERN-SPS and BNL-RHIC suggest 
what has been called sudden hadronization [|| or explo- 
sive behavior (jj|8|]- From the theoretical side, this phe- 
nomenon was recently associated with rapid changes in 
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the effective potential of QCD near the critical tempera- 
ture, such as predicted, for instance, by the Polyakov loop 
model , followed by spinodal decomposition || . 
Clearly, an understanding of the interplay between the 
typical space and time scales of the expanding plasma is 
welcome. Some att emp ts in this direction can be found 

in Refs. US El IH El 111 III US III 13 MM- 

In this paper, we discuss the mechanism responsible 
for the onset of instabilities in an expanding plasma. As 
a phenomenological model to mimic the case of the QGP, 
we use a relativistic plasma of quarks coupled to a chiral 
field. Although we derive a phenomenological chiral hy- 
drodynamics from a variational principle, we do not focus 
on the numerical solution of the resulting hydrodynamic 
transport equations. Rather, we consider the role played 
by chiral waves and pressure waves in the generation of 
unstable modes. We show that mechanical instabilities 
set in earlier than one would expect from the analysis of 
the thermodynamic potential decoupled from hydrody- 
namical modes. Therefore, the spinodal lines are shifted 
from the values one obtains studying only the chiral de- 
grees of freedom, modifying the phase diagram. 

It has been argued 0, IH E| that the first-order 
transition line which starts at the T = axis of the 
baryon chemical potential (/i) vs. temperature (T) phase 
diagram for QCD should come to an end at a critical 
point (/Ue,Te). Beyond this endpoint, for < fi < \xe, 
the QCD transition should become a smooth crossover, 
where spinodal instabilities do not occur. Results from 
lattice simulations at finite /j, also point in this direction 
pil I25I]. Here, we consider two illustrative scenarios for 
the temperature (T) vs. baryon chemical potential 
phase diagram in which a first order line starts at a point 
(/i,T) = (/x c ,0) and ends in a critical point, (^e,Te), 
with \jle < Mci beyond which, for < [i < [Ie, the phase 
transition becomes a smooth crossover. 

To model the mechanism of chiral symmetry break- 
ing present in QCD, we adopt a simple low-energy effec- 
tive chiral model: the linear cr-model coupled to quarks 
|26|. which in turn comprise the hydrodynamic degrees 
of freedom of the system. Similar approaches, relying on 
low-energy effective models for QCD and making use of 
a number of techniques to treat the expanding plasma, 
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can be found in the literature [3 13 13 13 13 The 
gas of quarks provides a thermal bath in which the long- 
wavelength modes of the chiral field evolve. The lat- 
ter plays the role of an order parameter in a Landau- 
Ginzburg approach to the description of the chiral phase 
transition 

The paper is organized as follows. In section II, we 
present our effective field theory model. In section III, 
we derive what we call phenomenological chiral hydro- 
dynamics from a variational formulation of relativistic 
hydrodynamics. In section IV, we study the onset of in- 
stabilities. In section V, we discuss our results. 



II. THE EFFECTIVE MODEL 

Let us consider a chiral field (f> = {a, it), where a is a 
scalar field and are pseudoscalar fields playing the role 
of the pions, coupled to two flavors of quarks according 
to the Lagrangian: 

C = q[h^ + - W(<f>)]q + \d^4> ~ V(<f>) . (1) 

Here q is the constituent-quark field q = (u,d) and 
fi q = /ii/3 is the quark chemical potential. The inter- 
action between the quarks and the chiral field is given 
by 

W(<f>) =g(a + ij 5 T-Tr) , (2) 

and 

v(4>) = ^(a 2 + n 2 -v 2 ) 2 -h q a (3) 

is the self-interaction potential for <f>. The parameters 
above are chosen such that chiral SUl(2) ® SUr{2) sym- 
metry is spontaneously broken in the vacuum. The vac- 
uum expectation values of the condensates are (a) = f n 
and (tt) = 0, where f v — 93 MeV is the pion decay con- 
stant. The explicit symmetry breaking term is due to 
the finite current-quark masses and is determined by the 
PC AC relation, giving h q = f^m^, where = 138 MeV 
is the pion mass. This yields v 2 = f 2 — m 2 /X 2 . The value 
of A 2 = 20 leads to a cr-mass, m 2 = 2A 2 / 2 + m 2 ., equal to 
600 MeV. In mean field theory, the purely bosonic part 
of this Lagrangian exhibits a second-order phase transi- 
tion |27j at T c = \/2v if the explicit symmetry breaking 
term, h q , is dropped. For h q ^ 0, the transition becomes 
a smooth crossover from the restored to broken symme- 
try phases. For g > 0, the finite-temperature one-loop 
effective potential also includes a contribution from the 
quark fermionic determinant |l3l Il4| . 

In what follows, we treat the gas of quarks as a heat 
bath for the chiral field, with temperature T and baryon- 
chemical potential /z. Then, one can integrate over the 
fermionic degrees of freedom, obtaining an effective the- 
ory for the chiral field <fi. To compute thermodynamic 



quantities, one needs the partition function 

Z = J[D(/)][Dq][Dq}e-^ /TdT ^ d3xCE , (4) 

where Ce is the Euclidean Lagrangian and V is the (infi- 
nite) volume of the plasma. Integrating over the fermions 
and using a classical approximation for the chiral field, 
we can write the thermodynamic potential as 

Q(T,n, <f) = V(cf>) - ^lndeti^ 1 + W{4>)]/T} , (5) 

where Ge is the fermionic Euclidean propagator satisfy- 
ing 

[fd T + if • V - ^7°] Ge(t,x;t',x') = 5 {i) {x~x') . 

(6) 

From the thermodynamic potential, (JSJ), one can ob- 
tain all the thermodynamic quantities of interest. The 
fermionic determinant that results from the functional 
integration over the quark fields can be calculated to one- 
loop order in the standard fashion |3 US • The effective 
Lagrangian for the chiral field, 0, in the presence of the 
quark thermal bath is then 

4// = 5^W-n(T,/i^). (7) 

III. CHIRAL HYDRODYNAMICS 

Given the thermodynamic potential, one can derive the 
total pressure and energy density, and obtain the con- 
served energy-momentum tensor, 7"^", for an expanding 
quark perfect fluid by using standard methods of ther- 
mal field theory (see, for instance, Ref. 01)- We prefer, 
instead, adopting an alternative approach to obtain the 
hydrodynamic equations for our system: the variational 
formulation |29l y0| • This approach provides a natural 
and unified way of merging chiral and fluid dynamics 
once the action of the system is specified. For a different 
treatment of the hydrodynamics of nuclear matter in the 
chiral limit, see [33j. 

The variational formulation of relativistic hydrody- 
namics has been recently applied to several physical sys- 
tems, especially in the realms of astrophysics and con- 
densed matter. In the former it has been used, for in- 
stance, to incorporate the effects of local turbulent mo- 
tion in supernova explosion mechanisms, and can prove 
to be useful in the analysis of the relativistic motion of 
blast waves in models for gamma-ray bursts |3l| . In soft 
condensed matter, it has been used to study the bubble 
dynamics in sonoluminescence experiments, in particular 
in deriving the relativistic generalization of the Rayleigh- 
Plesset equation |3|. In this section, we briefly review 
how to obtain the usual equations of relativistic hydro- 
dynamics within this framework, and derive in detail the 
case in which we have a chiral field coupled to an expand- 
ing quark fluid. 
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We describe the state of the fluid in terms of the four- 
velocity u^{x), the proper baryon density, n(x), and the 
proper entropy density, s(x). In what follows we assume 
that n and s are conserved, i.e.: 



d u (nu») = 
dp(au») = 



(8) 



The action that yields the hydrodynamic equations is 
given by [HOH 



(9) 



where e(n, s) is the proper energy density, from which 
temperature and chemical potential are obtained via the 
usual thermodynamic relations, T = de{n,s)/ds and 
/i = de(n,s)/dn. The equation for the hydrodynamic 
motion of the fluid is obtained by imposing the varia- 
tional principle with respect to n, s and u^, under the 
constraints JSJ) and the normalization condition 



1 



(10) 



These constraints can be incorporated in the variational 
principle in terms of Lagrangian multipliers, X(x), C( x ) 
and w(x), so that one imposes: 



S J d 4 x[ - e(n,s) + Xd l ,(nu f ") + Cd^su f ") 



2™ J ) 



= 0, 



(11) 



Equivalently, the fluid dynamics is given by the effective 
Lagrangian, £^ z / tlid ) = jrU l p d ) Sj ^ \, (, w), 

%( fluid) 



c 



eff 



-e(n,s) 



nu^d^X 
1 



- flti"0 M C+ -w(b\-1) . (12) 

Now the variables n, s, it M , A, £, and w are independent. 
Variations with respect to n, s and u 11 yield 



u^d^X = 



-T-vPdpC = 



-n<9 M A - sd u ( 







(13) 



while variations with respect to A, £ and w give simply 
the constraints JSJl and (fTTTfl - From (fTUf> and i fHjl . one 

obtains 



w = ri/U + Ts 
= e + P , 



(14) 



where p is the pressure, and identifies w as the enthalpy 
density. Using the Gibbs-Duhem relation, dp = sdT + 
ndn, and the property u 1/ (i9 AI u ly ) = 0, it is straightforward 
to obtain 

d v (wu^Uy) = nd a \i + sd^T + (d^u^) £u v 
= d u p + u v (d^Uv) £ 
= d uP , (15) 



which can be rewritten in the standard form 

d v T^ = 0, (16) 

where 



T^ v = wu^Uv - g^p 

= (p + e) u^u v - g^p 



(17) 



is the usual energy-momentum tensor of the fluid. 

It is important to notice that the effective Lagrangian, 
(|12f) . evaluated in the proper comoving frame of the fluid, 



(fluid) 
eff 



-e(n, s) + fj,n + Ts = p , 



(18) 



is nothing but minus the thermodynamic potential of 
the system. This fact will be useful to couple the chiral 
degrees of freedom to the quark fluid motion. 

The procedure outlined above for the derivation of rel- 
ativistic hydrodynamics presents a number of nice fea- 
tures. First, it can be easily generalized to include other 
degrees of freedom such as the chiral field. Secondly, 
although we have implicitly used a Minkowski metric, 
the generalization to the case of curved space-times is 
straightforward. This will be discussed later, when we 
consider applications to cosmology. Furthermore, from a 
practical point of view, once the variational approach is 
established one can use this method to obtain the optimal 
parameters for a given family of trial solutions. 

Let us now describe how to couple a chiral field to the 
fluid within the framework of the variational principle. 
The effective Lagrangian for the chiral field, Q), should 
be interpreted as the proper value of the Lagrangian 
where the thermalized quark fluid is at rest. There- 
fore, Eq. (|T%|) suggests that the part corresponding to 
the quark pressure should be replaced by expression 1|12|) 
in order to reconstruct the fluid motion of thermalized 
quarks. Apart from the constraints, we have the follow- 
ing action for the coupled system of the chiral field, </>, 
and the thermalized quark fluid motion: 

S tf+fiuid) = />a;J±(0 0(0/Y)_ e ( n>J}) 0)\ . (19) 



Here, the energy density, e, is related to the thermody- 
namic potential, f2, through the Legendre transformation 
e = Q + Ts + \m. Since 17 = — p, this is nothing but the 
general thermodynamic relation e = — p + Ts + fin. No- 
tice that now the energy density depends on the field 
variable (f> through Q = CI (T, fi, <j)) and, correspondingly, 
the second thermodynamic law should read 

de = Tds + fidn + Rdcj) . (20) 

The corresponding Gibbs-Duhcm relation becomes 

dp = sdT + ndfi - Rdcj) , (21) 

where we have defined a new quantity R (with four com- 
ponents Ri — d£l/dcf>i): 



R 



de(n, s, c 



4 



\ deb 



(22) 



The second line corresponds to Maxwell's relation and 
comes from Eq. Q21[l. R can be written as 



where 



d 3 k 1/E k (<f>) 



(2tt) 3 e [E k (4>)-n q ]/T + 1 



(23) 



-N) (24) 



is a scalar density for quarks. Here v q — 12 stands for 
the color-spin-isospin degeneracy of the quarks, Ek(4>) = 
{k 2 +m 2 q {^)) 1 / 2 , and m q {4>) = [g 2 cb 2 ) 1 / 2 = g{a 2 + tt 2 ) 1 / 2 
plays the role of an effective mass for the quarks. 

Thus, our coupled system is described by an effective 
Lagrangian 



C 



eff 



+fluU) = h d<j>) ( d ,><l>)-e(n,8,<l>)-nu> t 



su^d^C, + -w - 1) 



(25) 



The variation with respect to cb gives the equation of 
motion for the chiral field, 



□o 



-I? 



(26) 



The variation procedures with respect to the fluid vari- 
ables are almost exactly the same as before, except for 
the use of the new Gibbs-Duhem relation, 121|) . Thus, 
Eq. ((ToT) is modified to 



d v (wu^Uv) = d^p + Rd^cj), 
or, in the standard form, 

d v T vll = Rd^ , 



(27) 



(28) 



where T^ v is given by Eq. I|17|l and is not the total 
energy momentum tensor of the system, T^ v — — 
\d a <\>d a <\>g lLV + d^(j)d v (j) : which is of course conserved. 

Note that Eq. (|28() contains terms from the pure con- 
tribution of the field <j), such as V (4>) and dV/ dip, both 
in T^ v and in R. However, as we can see from the equiv- 
alent equation, l|27ll . some of these terms cancel out. To 
see more clearly the coupling between the fluid equation 
and the chiral field cb, it may be convenient to identify 
the fluid contribution of the pressure as 

pifiuid) = T i ndet [(G B 1 + W{<l>))/T] , (29) 



and the energy-momentum tensor of the fluid, 

= wu^ - v Ulmd) 9»v (30) 
Then the hydrodynamical motion is written as 

dt > T <jluM) _ pWmqdumqW = o . (31) 



IV. STABILITY ANALYSIS 

A numerical analysis of the system of differential equa- 
tions consisting of Q31(l and the equations of motion for 
a and tt was performed, for instance, in Refs. [TH Il4|. 
Here, we refrain from this approach and, instead, con- 
sider the behavior of small oscillations around equilib- 
rium to study the onset of instabilities. 

In what follows, we expand around equilibrium keeping 
terms up to first order in the perturbation, i.e.: 



n(x) = n eq + ni{x) , 

s(x) = s eq + si(ai) , 

<j){x) = <f> eq + ^(x) , 

u v (x) = u"+ul(x) . 



(32) 



Since the normalization of the four-velocity implies u\ 
(0, v\), and R eq = 0, the conservation laws (JSJ yield 



dni 
~~dt ' 

~dt 



n e „V ■ vi 



S eq V • Ul 



o, 

= , (33) 
and the chiral field equation of motion becomes 



where we have defined the (diagonal) chiral mass matrix 
~ d 2 e{n,s,cb) ] 

■ (36) 

The time component of Euler's equation, (|27|l . is au- 
tomatically satisfied given the conservation laws above, 
and its spatial components yield 



(36) 



For plane waves of the form ip\(x) — tpie lKx , where 
= (tu,k) and ipi stands for m, s\, <f>i or V\, we can 
rewrite the wave equation, (|34f> . and the spatial compo- 
nent of Euler's equation, (|36|l . as (dropping the tildes) 



(to 2 - k 2 

[u 2 - k 2 - ml) (Ji = —w eq R'vi 



m 2 ] 



0. 

k 



(uj 2 - p'k 2 ) V! = ukR'oi 

where we have defined 

' 'dp{e,n/s,(j>y 



P 



R' = 



dc 



dR a (c, n/s, 



de 



(37) 
(38) 
(39) 

(40) 

(41) 



eq 
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evaluated at constant (s/n) and 0, and used the conser- 
vation laws (|33[) in the form (dropping the tildes) 



ni 

n. 



(42) 



< <i 



One can see that pions decouple from the hydrody- 
namic modes, to first order, and play no major role in 
the following analysis. In this approximation, the inter- 
play between chiral modes and pressure modes happen 
entirely in the sigma sector of the chiral model. The 
dispersion relation then reads 



(u; 2 - p'k 2 )(io 2 -k 2 - ml) = w eq R' 2 k 2 



(43) 



For long wavelength fluctuations, we can approximate 
the roots of (g3J> by 



w ea R'' 



>>>n 



0{k 2 ) , 



0{k 2 ) 



(44) 



(45) 



These waves can be identified as pressure modes of 
frequency lu p , moving with the speed of sound 



w ea R L 



1/2 



and chiral modes with frequency Lu a and mass 



d 2 e(n, s, a) 



da 2 



(47) 



ft; 



The onset of instabilities takes place when uj 2 < 5] , 
i.e., either for v 2 ound < or for m 2 < 0. We see from 
Eq. (|46H that pressure modes become unstable before chi- 
ral modes do. 

To study the onset of instabilities within the phase dia- 
gram of the low-energy effective model for QCD described 
previously, we consider two illustrative cases. The first 
corresponds to a coupling between quarks and the chiral 
field small enough for the transition at /i = to be a 
smooth crossover. We use g = 3.3, which yields a con- 
stituent quark mass of 307 MeV, about 1/3 of the nucleon 
mass. In this case there is a line of first-order transitions 
starting at T = 0, \i c = 920 MeV, and ending at a criti- 
cal point E with T E = 98 MeV and \i E = 630 MeV. In 
the second scenario we take a larger coupling, g = 5.5, so 
that even at [i = there is a first-order phase transition, 
in this case at T r w 123 MeV [ll 111 111. 



V. DISCUSSION OF RESULTS 

Let us now discuss the onset of instabilities in the 
phase diagram of our effective model. In Fig. Qwe plot 
the phase diagram in the (//, T)-plane for the couplings 
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FIG. 1: Phase diagram for the effective model in the (fj,, T)- 
plane, for g — 3.3. The coexistence line ending at the critical 
point E is represented by the full curve. The spinodal lines 
(mj; = 0) are shown as dashed curves. 
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FIG. 2: Phase diagram for the effective model in the (n, T) 
plane, for g — 3.3. Line conventions are the same as in Fig. 1. 
Inside the full curve the system is mechanically unstable. 



g = 3.3. The coexistence line is represented by the solid 
curve. The lines where pressure waves become unstable 
(v soun d — 0) are shown as dashed curves. They corre- 
spond to the supercooling and superheating spinodals. 

The phase diagram in the (n, T)-plane is shown in 
Fig. [21 The phase border of the coexistence region is 
represented by the solid line and dashed lines stand for 
the spinodals. The sector between the lines on the right 
of the critical point E corresponds to supercooled states 
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FIG. 3: The same as Fig. [T]for g = 5.5. 
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FIG. 4: The same as Fig. Hfor g = 5.5. 



in the chirally symmetric phase. Superheated states cor- 
respond to the area on the left. The domain inside the 
dashed lines corresponds to mechanically unstable states 
which undergo spinodal decomposition. 

Results for the second choice of the coupling constant, 
g = 5.5, are shown in Figs. and 0| As we have men- 
tioned, the difference of this case with respect to the first 
is that the coexistence line now reaches /x = 0. Because 
of this, the spinodals cross each other in the (n, T) plane, 
and part of the phase coexistence region can be occupied 
by both supercooled and superheated states. Below the 
crossing no metastable states exist. 

We have shown that the process of domain formation 
and growth in the phase conversion of strongly interact- 
ing matter from a chirally symmetric to a chirally bro- 
ken state in an expanding QGP is dominated by the ex- 
ponential increase of hydrodynamical sound-like modes, 
while chiral-like modes remain stable. We have mapped 
the boundaries of the unstable region for two illustrative 
cases, determining how much supercooling (or superheat- 
ing) is necessary for the onset of spinodal instabilities. 

In the case of relativistic heavy ion reactions, the short 
time scale of expansion and effects due to the finite size 
of the system could change significantly our results. In 
order to address these points one should perform exten- 
sive numerical simulations to study the development of 
instabilities in different scenarios. Results in this direc- 
tion will be presented elsewhere. 



Acknowledgments 



We thank A. Dumitru, K. Paech, D. Rischke, J. 
Schaffner-Bielich and H. Stocker for fruitful discussions. 
This work was partially supported by CAPES, CNPq, 
FAPERJ and FUJB/UFRJ. 



[1] F. Karsch, Nucl. Phys. A 698, 199 (2002). 
[2] J. Harris and B. Miiller, Ann. Rev. Nucl. Part. Sci. 46, 
71 (1996). 

[3] H.G. Ritter and X.N. Wang (ed), Proceedings of Quark 
Matter 2004, J. Phys. G 30 S633-S1425 (2004). 

[4] E.W. Kolb and M.S. Turner, The Early Universe 
(Addison- Wesley, Redwood City, 1990). 

[5] J.D. Gunton, M. San Miguel and P.S. Sahni in Phase 
Transitions and Critical Phenomena (Eds.: C. Domb and 
J. L. Lebowitz, Academic Press, London, 1983), v. 8. 

[6] J. Rafelski and J. Letessier, Phys. Rev. Lett. 85, 4695 
(2000). 

[7] A. Dumitru and R.D. Pisarski, Nucl. Phys. A 698, 444 



(2002). 

[8] O. Scavenius, A. Dumitru and A.D. Jackson, Phys. Rev. 

Lett. 87, 182302 (2001). 
[9] R.D. Pisarski, Phys. Rev. D 62, 111501 (2000); A. Du- 
mitru and R. D. Pisarski, Phys. Lett. B 504, 282 (2001). 

[10] L.P. Csernai and I.N. Mishustin, Phys. Rev. Lett. 74, 
5005 (1995); A. Abada and J. Aichelin, Phys. Rev. Lett. 
74, 3130 (1995); A. Abada and M. C. Birse, Phys. Rev. 
D 55, 6887 (1997); 

[11] I.N. Mishustin and O. Scavenius, Phys. Rev. Lett. 83, 
3134 (1999). 

[12] O. Scavenius and A. Dumitru, Phys. Rev. Lett. 83, 4697 
(1999). 



7 



[13] O. Scavenius, A. Dumitru, E.S. Fraga, J. T. Lenaghan 

and A. D. Jackson, Phys. Rev. D 63, 116003 (2001). 
[14] K. Paech, H. Stoecker and A. Dumitru, Phys. Rev. C 68, 

044907 (2003); K. Paech and A. Dumitru, Phys. Lett. B 

623, 200 (2005). 
[15] E.E. Zabrodin, L.V. Bravina, L.P. Csernai, H. Stocker 

and W. Greiner, Phys. Lett. B 423, 373 (1998). 
[16] P. Shukla and A. K. Mohanty, Phys. Rev. C 64, 054910 

(2001). 

[17] C. Spieles, H. Stocker and C. Greiner, Phys. Rev. C 57, 
908 (1998). 

[18] E.S. Fraga and R. Venugopalan, Physica A 345, 121 

(2005); Braz. J. Phys. 34, 315 (2004). 
[19] J. Randrup, Phys. Rev. Lett. 92, 122301 (2004); J. Ran- 

drup, nucl-th/0406031 V. Koch, A. Majumder and 

J. Randrup, nucl-th/0509030 
[20] E.S. Fraga and G. Krein, Phys. Lett. B 614, 181 (2005). 
[21] MA. Stephanov, K. Rajagopal and E.V. Shuryak, Phys. 

Rev. Lett. 81, 4816 (1998); Phys. Rev. D 60, 114028 

(1999). 

[22] M. A. Halasz, A. D. Jackson, R. E. Shrock, 
M. A. Stephanov and J. J. Verbaarschot, Phys. Rev. D 
58, 096007 (1998); J. Berges and K. Rajagopal, Nucl. 
Phys. B 538, 215 (1999); T.M. Schwarz, S.P. Klevansky 



and G. Papp, Phys. Rev. C 60, 055205 (1999). 
[23] O. Scavenius, A. Mocsy, I. N. Mishustin and 

D.H. Rischke, Phys. Rev. C 64, 045202 (2001). 
[24] Z. Fodor and S.D. Katz, Phys. Lett. B 534, 87 (2002); 

JHEP 0203, 014 (2002). 
[25] C.R. Allton et al, Phys. Rev. D 66, 074507 (2002). 
[26] M. Gell-Mann and M. Levy, Nuovo Cim. 16, 705 (1960); 

R. D. Pisarski, Phys. Rev. Lett. 76, 3084 (1996). 
[27] R.D. Pisarski and F. Wilczek, Phys. Rev. D 29, 338 

(1984). 

[28] J. Kapusta, Finite Temperature Field Theory (Cam- 
bridge University Press, Cambridge, 1989). 

[29] M. Bailyn, Phys. Rev. D 22, 267 (1980). 

[30] H.T. Elze, Y. Hama, T. Kodama, M. Makler and J. Rafel- 
ski, J. Phys. G 25, 1935 (1999). 

[31] H. A. R. Goncalves, S. B. Duarte, T. Kodama and V. 
DAvila, Astron. & Space Science 194, 313 (1992); T. 
Kodama, R. Donangelo and M. Guidry, Int. J. Theor. 
Phys. C 9, 745 (1998). 

[32] H.T. Elze, T. Kodama, J. Rafelski, Phys. Rev. E 57, 4170 
(1998). 

[33] D.T. Son, Phys. Rev. Lett. 84, 3771 (2000). 



